Relaxation properties in a lattice gas model with asymmetrical particles 
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We study the relaxation process in a two-dimensional lattice gas model, where the interactions come 
from the excluded volume. In this model particles have three arms with an asymmetrical shape, 
which results in geometrical frustration that inhibits full packing. A dynamical crossover is found 
at the arm percolation of the particles, from a dynamical behavior characterized by a single step 
relaxation above the transition, to a two-step decay below it. Relaxation functions of the self-part 
of density fluctuations are well fitted by a stretched exponential form, with a j3 exponent decreasing 
when the temperature is lowered until the percolation transition is reached, and constant below it. 
The structural arrest of the model seems to happen only at the maximum density of the model, 
where both the inverse diffusivity and the relaxation time of density fluctuations diverge with a 
power law. The dynamical non linear susceptibility, defined as the fluctuations of the self-overlap 
autocorrelation, exhibits a peak at some characteristic time, which seems to diverge at the maximum 
density as well. 



I. INTRODUCTION 

Most glassy systems such as structural glasses, ionic 
conductor, supercooled liquids, polymer, colloid, and 
spin glasses jjj have similar complex dynamical behav- 
ior. As the temperature is lowered the relaxation times 
increase drastically, and the relaxation functions deviate 
strongly from a single exponential function at some tem- 
perature T* well above the dynamical transition. In the 
long time regime they can be well fitted by a stretched 
exponential or Kohlrausch- Williams- Watts Q function 
/(*) = f exp [-{t/rf) , with < < 1. 

There are two mechanism driving to non exponential 
relaxation. In disordered model like spin glasses, it is 
caused by the existence of unfrustrated ferromagnetic- 
type clusters of interactions and therefore is a direct 
consequence of the quenched disorder || . Another mech- 
anism in frustrated systems is based on the percolation 
transition of the Kasteleyn-Fortuin and Coniglio-Klein 
cluster Q. In this approach disorder is not needed to 
obtain nonexponential relaxation 

In spin glasses there is a thermodynamic transition at 
a defined temperature TsG^ where the non linear suscep- 
tibility diverges. As T$g is approached, the static corre- 
lation length becomes larger. In glass forming liquids it 
seems that there is no sharp thermodynamic transition, 
and no diverging static length. However numerical stud- 
ies have identified long lived dynamical structures which 
are characterized by a typical length and a typical relax- 
ation time, which depend on temperature and density || . 
In order to characterize this behavior a dynamical non- 
linear susceptibility was introduced by Donati et al. [jlO| , 
both for spin models and for structural glasses. They 



have found that the dynamical susceptibility exhibits a 
maximum at some characteristic time which diverges as 
the dynamical transition temperature is approached from 
above. Also in the annealed version of the frustrated 



Ising lattice gas model 11 1 a similar behavior has been 
found, while in the quenched version the dynamical sus- 
ceptibility is always increasing, due to the presence of the 
thermodynamic transition |l2| . 

The understanding of the macroscopic process of relax- 
ation in structural glasses, starting from the microscopic 
motion of the particles, has been boarded using differ- 
ent microscopic models, as for example the hard square 
model pq ], the kinetically constrained models Jl4|,15|, 
and the frustrated Ising lattice gas (FILG) model [11|. 
These models have reproduced some aspects of the glassy 
phenomenology, and recently the FILG model has been 
studied in its quenched and annealed version with kinetic 
constraints Jl2]. The results found in the annealed ver- 
sion seems to be closer to the experimental ones. In this 
paper we consider a two-dimensional grometrical model, 
which contains as main ingredients only geometrical frus- 
tration without quenched disorder and without kinetic 
constraints, as quenched disorder is not appropriate to 
study structural glasses and kinetic constraints are some 
how artificial. Similar models have already been pro- 
posed jll| and applied to study granular material (l(| . 

In Sec. II we present the model and in Sec. Ill we 
study its percolation properties. In Sec. IV we show the 
dynamical results, and finally in Sec. V we present our 
conclusions. 
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II. THE MODEL 

In this paper we introduce a model which can be con- 
sidered as an illustration of the concept of frustration 
arising as a packing problem. In systems without un- 
derlying crystalline order, frustration is typically gener- 
ated by the geometrical shape of the molecules, which 
prevents the formation of close-packed configurations at 
low temperature or high density; for systems with un- 
derlying crystalline order, frustration arises when the lo- 
cal arrangement of molecules kinetically prevents all the 
molecules from reaching the crystalline state. 

An example of glass-former that has difficulty in 
achieving crystalline order is the ortho-terphenyl, whose 
molecule is made of three rings. This system can be 
loosely modeled with a simple lattice model, in which 
"T" shaped objects occupy the vertices of a square lat- 
tice with one of four possible orientations. 

Assuming that the arms cannot overlap due to ex- 
cluded volume, we see that only for some relative ori- 
entations two particles can occupy nearest-neighbor ver- 
tices. Consequently, depending on the local arrangement 
of particles, there are sites on the lattice that cannot be 
occupied (see Fig. Q). This type of "packing" frustration 
thus induces defects or holes in the system. This model 
resembles the hard-square lattice gas model [|l3| which 
can be seen as "+" shaped objects on the vertices of a 
square lattice with excluded volume interaction. A very 
important difference between these two models is the in- 
ternal degree of freedom due to the particle shape which 
is absent in the latter. 

We consider a two dimensional square lattice and im- 
pose periodic boundary conditions. In our system the 
maximum of density is p max = 2/3 at which all possible 
bonds are occupied by an arm. A configuration of density 
Pmax is a ground state of the system, corresponding to 
chemical potential /i-nxior temperature T — ► 0. It can 
be obtained for any size constructing larger systems from 
smaller ones with /9 max and appropriated boundary con- 
ditions. In this way one can build an extensive number 
of different ground states that lack spatial order. 

We have simulated the diffusion and rotation dynamics 
of this model by Monte Carlo methods. The dynamics 
of the particles is given by the following algorithm: i) 
Pick up a particle at random; ii) Pick up a site at ran- 
dom between the four nearest neighbor ones; iii) Choose 
randomly an orientation of the particle; iv) If it does not 
cause the overlapping of two arms, move the particle in 
the given site with the given orientation; v) If the diffu- 
sion movement is not possible, choose a random orienta- 
tion and try to rotate the particle to this new orientation; 
vii) Advance the clock by 1/N S , where N s is the number 
of sites, and go to i). 

Studying the dynamics by Monte Carlo simulations, 
the finite size effects are larger when we are near to p max , 



because the particles can be enclosed in cages, and the 
diffusion is blocked. In lattices of finite size, cages may 
be indefinitely stable. This has also been observed in the 
hard-square lattice gas model ]l3] | where at the maximum 
of density the particles are on the diagonals, but in lat- 
tices of finite size cages of particles indefinitely stable are 
formed at lower densities. 



III. PERCOLATION TRANSITION 

In order to investigate whether the percolation transi- 
tion has effects on the dynamics, in this section we ana- 
lyze the percolation properties of the model, and relate 
the percolation density with a change in the dynamical 
properties of the model. The particles have three arms, 
and there are two bonds per site on the lattice, so the den- 
sity of bonds occupied by an arm is given by a = Sp/2, 
where p is the density of particles. Therefore, if the corre- 
lations of the arms were not important, the arm percola- 
tion would occur at the density a c = 1/2, corresponding 
to p c = 1/3. Nevertheless we expect some correlation 
effects. 

We have simulated our system for various lattice sizes 
around the percolation density in order to determine p c . 
For each density we have reached equilibrium and then, 
taking 10 4 steps, we have evaluated the probability of 
existence of a spanning cluster P and the mean cluster 
size S = s 2 n s , where n s is the density of finite clusters 
of size s. 

Around the percolation density the averaged quantities 
P(p) and S(p), for different values of the lattice size L, 
should obey to the finite size scaling p8| 



P(p)=F P \l}l»{p-p c ) 



s(p)^l->^f s \l 1 ^(p- Pc ) 



(1) 

(2) 



where 7 and v are then critical exponents, and Fp(x) and 
Fs(x) are universal functions of an adimensional quantity 
x. Fig. U shows the finite size scaling of P and S. We 
have selected the values of the exponent corresponding 
to the two-dimensional site-bond percolation universality 
class psf , that is v = 4/3 and 7 = 43/18 and we have 
found p c = 0.315 ± 0.003. The critical density can also 
be determined from P(p) as the density at which curves 
corresponding to different sizes cross. From the inset of 
Fig. ||(a) we see that p c is between 0.315 and 0.32, so that 
the arm correlations and the thermal process decrease the 
critical density with respect to the random bond problem. 



IV. DYNAMICAL RESULTS 

We wish now to define a self-overlap parameter, in or- 
der to measure the autocorrelation of the model. Dif- 
ferent definitions of q can be used, and the form of the 
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relaxation functions is different for different overlap pa- 
rameters. Here we define a self-overlap parameter which 
takes into account the orientation of the particle, rather 
that its position. The orientation of the particle is de- 
fined by the discrete values of the angle, fa = 0, n/2, n, 
or 3n/2, and we define the self-overlap 

q(t) = — J2ni(t')ni(t' + t) cos [<i>i{t' + t)-<Pi(t')] (3) 

i 

where rii(t) = 0, 1 is the occupation number of site i at 
time t, fa (t) is the orientation of the particle on site i at 
time t, and N is the number of sites. This parameter is 
a generalization of the self-overlap defined in Ref. p2| , 
where the orientation fa plays the role of the spin vari- 
ables. When all particles are frozen then q = 1. The as- 
sociated dynamical nonlinear susceptibility x(i) is given 
by 

X (t) = N[(q(t) 2 )-( q (t)) 2 ], (4) 

where the average (• • ■) is done on the reference time t' . 

In Fig. H we show the relaxation function of the self- 
overlap parameter, for a system of size 64 2 and densities 
between p = 0.2 and 0.66. Each curve is obtained aver- 
aging over a time interval for t' of 10 4 — 10 8 Monte Carlo 
steps, depending of the density. For high and intermedi- 
ate densities we observe a two step relaxation functions, 
which gives evidence of two well separate time scales in 
the system. The first short time decay of the relaxation 
functions is due to the rotations of the particles in a 
frozen environment, which appears as quenched on this 
time scale, while the second decay is due to the evolution 
of the environment and final relaxation to equilibrium. 

For very low densities p < 0.1 the relaxation function 
is well fitted by an exponential form. From p » 0.1 to 
p « 0.32 we can fit the relaxation function by a stretched 
exponential form, with the exponent (3 decreasing slowly 
until j3 ~ 0.9 at p « 0.32. For higher values of densities 
it is impossible to fit the relaxation function with a sin- 
gle form over the whole time interval. As we have seen, 
p w 0.32 corresponds to arm percolation. This suggests 
that the percolation transition induces a crossover in the 
dynamics, from a single relaxation regime to two tempo- 
ral scales for the relaxation of the particles. At densities 
higher than p w 0.32, the long time tail of the relax- 
ation functions is well fitted by a stretched exponential 
form, where the exponent (3 depends very weakly from 
the density (it is constant within the errors) and ranges 
between (3 = 0.64 and (3 = 0.71. We show in Fig. | 
the time-density superposition of the relaxation function 
(q(t)), for densities between p — 0.62 and p — 0.66. Also 
shown as a dashed line is a fit with a stretched exponen- 
tial function. 

In Fig. we show the intermediate time behavior of 
the self-overlap relaxation functions, for densities p — 



0.62, 0.64 and 0.65. We have tried to fit them with the 
simplified form of the function predicted by the MCT, 

(q(t)) = / + At a - Bt b (5) 

where the fitting parameters are /, A, B and A, while a 
and b are related to A by the transcendental equation 

r 2 (i-«) = r 2 (i + 6) 

r(i-2o) r(i + 26) ■ {> 

The values of A given from the fits are constant within 
the errors, with a mean value A = 0.785 ± 0.005, which 
correspond to exponents a — 0.285±0.005 and b — 0.50± 
0.01. 

In Fig. |^ we show the dynamical non linear suscepti- 
bility for some values of the density. The maximum in 
the susceptibility x(i*) and the time of the maximum t* 
seem to diverge together when the density grows. This 
has also been found previously in other models as p-spin, 
Lennard- Jones binary mixture fTpf , and in the annealed 
version of the frustrated Ising lattice gas model . In 
our model we obtain that the maximum of x(^*) can be 
fitted by the power law x(i*) oc (/5 max — p)~ a - Here we 
have p max = 0.664 ± 0.002 and a = 0.71 ± 0.02. The 
equilibrium value is x(t ~ * 00 ) — P 2 /2 for low densities 
and x(t — > oo) = 1/2 for the higher ones. 

The density-density autocorrelation function and its 
dependence with the time is an important property which 
characterizes the glassy behavior. We have studied the 
self-part of the autocorrelation function of the density 
fluctuations defined as 

^w = ^E elk(ri(t ' +t) " r,(t ' )) )' ( ? ) 

where Ti(t) is the position of the ith particle in units of 
the lattice constant. The wave vector can take the values 
k = (27r/£)n, where n has integer components n x and 
n y ranging from to L/2. 

Fig. shows F£(t) corresponding to k x — tt and k y = 
for different densities. For all densities the whole time 
interval of the autocorrelation function can be fitted by 
a stretched exponential function, f(t) — exp [— (t/r)' 3 ], 
where the exponent (3 depends of the density. In Fig. 
H we show j3 as a function of the density. We can see 
that the exponent (3 decreases with the density until a 
density near p « 0.32 is reached. From this density the 
exponent becomes constant (3 « 0.82 (within of the error 
bars). At densities near to p max the finite size effects be- 
come important and (3 deviates from the constant value. 
As in the a-relaxation of the self-overlap parameter, the 
effect of the percolation on the relaxation of F£(t) is to 
leave (3 constant, but now there is not two step relaxation 
process. 

If we look at the relaxation of the position of particles 
there is no vibration relaxation, the relaxation begins 
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when frozen domains evolve. Instead the relaxation of 
the orientation of the particles has a relaxation process 
where it is possible the vibration of the rotational degree 
of freedom. Nevertheless after the arm percolation the 
position and orientation relaxation in the a-relaxation 
have an stretched exponential form with the exponent (3 
constant (within the errors) although with different value 
each other. 

The relaxation time r is obtained from the fitting of 
F£(i) with a stretched exponential function, and can be 
fitted by a power law r cx (p max — p) -7 , with 7 = 2.70 ± 
0.02 and p max = 0.666 ± 0.003 (see Fig. |). Note that 
the value of 7 coincides within the errors with the one 
deduced by the exponents a and b given by the fit of Fig. 
H using the relation predicted by Mode-Coupling Theory 



which gives 7 = 2.74 ± 0.03. 

We have calculated the diffusion coefficient from the 
mean-square displacement (Ar(t) 2 ) at very long times. 
The values obtained for D are well fitted by a power law 
near of p max , D cx (p max — p) 7 with 7 = 2.70 ± 0.02 and 
Pmax = 0.666 ±0.003 (see the inset of Fig. §). The power 
law agrees with the law found in some lattice gas models 
(the kinetically constrained model but it is different 
of the law found in the hard-square lattice gas model 
jn|. This singular behavior of D is in accordance with 
the prediction of MCT. Form the behavior of D and t wc 
arrive to D^ 1 cx r and so to the Stokes- Einstein relation 
D^ 1 cx r\ (due to t cx 77). 

V. CONCLUSIONS 

We have proposed a two-dimensional geometrical mo- 
del, based on the concept of geometrical frustration which 
is generated by the particle shape. This model has nei- 
ther quenched disorder nor kinetic constraints. Percola- 
tion in the model has been studied, and it has been found 
that at the critical density of arm percolation the model 
shows a dynamical crossover, characterized by the onset 
of a two step relaxation in the orientational relaxation of 
the particles. Below the transition, the long time regime 
of orientational and positional relaxation functions can 
be fitted by a stretched exponential with an exponent (3 
constant within the errors. The dynamical results have 
been compared with the prediction of the MCT, find- 
ing good agreement. A self-overlap parameter has been 
defined which takes into account the orientation of par- 
ticles, and the corresponding dynamical non linear sus- 
ceptibility has been studied. It has been found that this 
dynamical susceptibility shows a peak at some character- 
istic time, which gives evidence of long lived dynamical 
structures with a growing length and relaxation time, as 
found in some molecular dynamics simulations. 
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FIG. 1. Schematic picture of one particular configuration 
in a system size of 4 2 and density p = 9/16. 
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FIG. 2. Finite size scaling of (a) P(p) and (b) S(p) for 
lattice sizes L = 40, 50, 60, 70, 80, 90, 100, 110, and 120. 
Inset: (a) P(p) and (b) S(p) for the same sizes. 
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FIG. 3. Relaxation functions for the self-overlap for system 
size 64 2 and densities p = 0.2, 0.3, 0.4, 0.5, 0.55, 0.6, 0.62, 
0.63, 0.64, 0.65, 0.655, 0.657, and 0.66. 
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FIG. 4. Time-density superposition principle for the relax- 
ation function of the self-overlap, for densities p = 0.62, 0.63, 
0.64, 0.65, 0.655, 0.657, and 0.66. The dashed line is a fitting 
function corresponding to a stretched exponential form with 
exponent /3 = 0.71. 
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FIG. 5. Fit of the intermediate time regime of the relax- 
ation functions of the self-overlap, for densities p = 0.62, 0.64 
and 0.65. The fitting function is / + At~ a - Bt b (dotted line), 
where the fitting parameters are /, A, B, and A, and a and b 
are given by the relation (6) . 



FIG. 6. Dynamical susceptibility for L — 64 and densities 
p = 0.6, 0.62, 0.63, 0.64, 0.65, 0.655, 0.657, and 0.66. In- 
set: the maximum x(t*) as a function of density. The fitting 
function is a power law x(0 = 0.12(0.664 — p)~ 0,7 . 
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FIG. 7. Relaxation functions for the self-part of the den- 
sity-density autocorrelation function for k x — ir and k y = 
and the same system size and densities of Fig. 2. Dotted lines 
are fitting functions corresponding to stretched exponential 
functions. 
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FIG. 8. Parameter j3 as a function of the density, obtained 
by fitting F£(t) for k x — tv and k y — with the function 
f(t) = exp [— (t/r) 13 ] . The densities represented are p = 0.05, 
0.1, 0.15, 0.2, 0.25, 0.3, 0.33, 0.35, 0.37, 0.4, 0.45, 0.5, 0.55, 
0.6, 0.63, 0.64, 0.65, 0.655, and 0.66. The solid line corre- 
sponds to (3 = 0.82. 




FIG. 9. The relaxation time T a as a function of density p. 
Solid circles: relaxation times obtained from the fit in Fig. 7 
Asterisks: relaxation times obtained from the fit in Fig. i 
The solid lines are power law functions r a = a(0.666 — p)~ 2 ' , 
with a — 0.115 in the former case and a — 0.04 in the latter. 
Inset: the diffusion constant D as function of density p. The 
fitting function is a power law D = 8(0.666 — p) 2 ' 7 . 
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